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Abstract. A Feshbach resonance in the s-wave scattering length occurs if the energy 
of the two atoms in the incoming open channel is close to the energy of a bound state 
in a coupled closed channel. Starting from the microscopic hamiltonian that describes 
this situation, we derive the effective atom-molecule theory for a Bose gas near a 
Feshbach resonance. In order to take into account all two-body processes, we have 
to dress the bare couplings of the atom-molecule model with ladder diagrams. This 
results in a quantum field theory that exactly reproduces the scattering amplitude of 
the atoms and the bound-state energy of the molecules. Since these properties are 
incorporated at the quantum level, the theory can be applied both above and below 
the critical temperature of the gas. Moreover, making use of the true interatomic 
potentials ensures that no divergences are encountered at any stage of the calculation. 
We also present the mean-field theory for the Bose-Einstein condensed phase of the 
gas. 
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1. Introduction 

The recent experiment by Donley et al. [|T| has made it clear that near a Feshbach 
resonance the coherence between atoms and molecules can have a profound role in the 
dynamics of an atomic Bose-Einstein condensate. As a result, it has become an urgent 
problem to understand from first principles how to properly incorporate the possible 
coherence between atoms and molecules into the theory of an interacting Bose gas. 
In the last four years important progress towards a solution of this problem has been 
made by a number of groups P, [3|, |], ||, ^| [7]]. Nevertheless, it appears that a fully 
satisfactory theory, which obeys all the requirements that on physical grounds can a 
priori be imposed upon the theory, still needs to be developed. It is with this goal in 
mind that the present contribution has been written. 

A first requirement for the theory is that it is based on an adequate microscopic 
description of a Feshbach resonance. Feshbach resonances are already known for a long 
time in nuclear physics ||, but have only more recently been predicted to occur in 
ultracold atomic gases || [L(| . The defining feature of such resonances is that they can 
only occur in a multi-channel scattering problem. More precisely, a Feshbach resonance 
occurs when the kinetic energy of the particles in the incoming open channel is equal to 
the energy of a bound state in a closed channel that is (weakly) coupled to the incoming 
channel. It is important for our purposes that the physics of a Feshbach resonance is 
quite different from the physics of a resonance in a single-channel scattering problem 
with a bound state near the continuum threshold of the particles. This can most easily 
be seen from the fact that in the latter case the bound state in the incoming channel is 
only very weakly bound and the extent of its wave function is, therefore, always much 
larger than the range of the interaction between the particles. For a Feshbach resonance 
the extent of the bound state wave function is generally of the same size as the range of 
the interactions, because it usually corresponds to a deeply bound state in the potential 
of another and closed channel. As a result also the quantum numbers of the bound state 
are different from the quantum numbers of the incoming particles, which is clearly not 
the case for the single-channel scattering problem. 

A second requirement follows from the fact that the many-body theory based on the 
microscopic description envisaged above is free of ultra-violet divergences at any level 
of approximation, and in particular, for both the normal and superfluid phases of the 
atomic Bose gas. The reason for this is that microscopically the interatomic interactions 
responsible for the Feshbach resonance are determined by short-range potentials, which 
cut off all the momentum integrals that arise when we diagramatically include the effects 
of the interactions. As a result the desired effective atom-molecule hamiltonian of the 
gas can contain only terms with coefficients that are finite and cut-off independent. 
Moreover, the coefficients are determined by only a small number of experimental 
parameters such as the position and the width of the Feshbach resonance, for instance. 

The most crucial additional requirement on the effective atom-molecule hamiltonian 
is that it must exactly reproduce the two-body physics of the atomic gas, i.e., by solving 
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the hamiltonian in the Hilbert space of two atoms we must recover the correct scattering 
amplitude of the atoms and also the correct binding energy and quantum numbers of the 
molecule. Note that in this manner the molecular properties are exactly incorporated at 
the quantum level, which is important for several reasons. First, the theory can now be 
applied both in the normal and in the superfluid phase of the gas. Second, we are not 
restricted to a mean-field description of the gas and are now also able to systematically 
study fluctuation effects, which we have recently shown to be of importance under 
the experimental conditions of interest ||11|| . Finally, the theory can be immediately 
generalized to atomic Fermi gases near a Feshbach resonance, which are of great current 
interest in view of the prospect of creating new neutral BCS superfluids |13|, 
Also in this case fluctuation effects are known to be essential |15|, [16| and can only be 
accounted for after the molecular properties are exactly incorporated in the quantum 
theory. 

As already implicitly mentioned, the effective atom-molecule hamiltonian is in first 
instance the most convenient way to arrive at a mean-field description of the gas. In 
the case of a Bose-Einstein condensed atomic gas, an exact property of the system is 
that it has a gapless excitation. Therefore, any physically reasonable mean-field theory 
must have the same property. From a fundamental point of view the gapless excitation 
is due to the fact that a Bose-Einstein condensate spontaneously breaks the global U{1) 
symmetry associated with the conservation of the total number of atoms. For the mean- 
field theory to fulfill this requirement automatically it must be formulated such that it 
does not contain an anomalous density or pairing field [|17|, [18| . Including an anomalous 
density in general leads to a double counting of the interaction effects and, therefore, 
destroys the gaplessness of the mean-field theory. 

In agreement with the above discussion, the paper is organized as follows. In 
section H we start from a microscopic hamiltonian for an atomic Bose gas with a Feshbach 
resonance and derive from this a bare atom-molecule theory. To make contact with 
the experimentally known parameters of the Feshbach resonance, we then carry out in 
section |3] a complete ladder summation to arrive at the desired effective atom-molecule 
hamiltonian that incorporates exactly all the relevant two-atom physics. In section |] we 
show in particular how the correct properties of the molecule are recovered. After that 
we discuss in section [5] how to arrive at the simplest mean-field theory that is appropriate 
at the low temperatures of interest experimentally, where the thermal component of the 
gas can be neglected. We finally end in section || with our conclusions. 



2. Bare atom-molecule theory 

Without loss of generality we can consider the simplest situation in which a Feshbach 
resonance arises, i.e., we consider a homogeneous gas of identical atoms with two internal 
states denoted by | |) and | [). The atoms in the two states interact via the potentials 
V||(x — x') and Vj^(x — x'), respectively. The state | {) has an energy AfiB/2 with 
respect to the state | |) due to the Zeeman interaction with the magnetic field B. The 
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coupling between the two states, which from the atomic physics point of view is due to 
the hyperfine interaction, is denoted by V^(x — x'). Putting everything together our 
microscopic hamiltonian is thus given by 



H 



dxip^ 



n 2 v 2 1 

2m + 2 



cixVI(x / )K TT (x-x / )^ T (x , ) 



^ T (x) 



dxip] (x) 



+ ~Y~ + n I rfxVj(x')ni(x " x')^(x') 



1 



+ 2 MX 



dk [^(x)^}(x')K n (x - x')^(x')^(x) + h.c. 



where the potential Vj_j_(x — x') is assumed to contain the bound state responsible for 
the Feshbach resonance. Using a Hubbard- Stratonovich transformation to decouple 
this part of the hamiltonian 20|, we introduce the molecular field operator ^ m (x) 
that annihilates a molecule at position x. In the approximation that we are close to 
resonance, only a single bound state contributes and this operator has the property that 

(^(x)^(x')) = >/2$ m ((x + x')/2)) Xm (x - x'). (2) 

The properly normalized and symmetrized bound state wave function in the potential 
Vjj,(x — x'), which we choose to be real for simplicity, obeys the Schrodinger equation 



ti 2 V 2 



m 



(3) 



After the Hubbard- Stratonovich transformation we obtain the bare hamiltonian for the 
coupled atom-molecule system. It reads 

h 2 v 2 l 



H 



dxipl 



+ / dx^(x) 



2m 

h 2 v 2 



dkM(x.')V n {x-k)ip^(k) 



^t(x) 



+ AfiB + E n 



+ / dx / dx 



(4) 



4m 

g n (x - x')A((x + x')/2)^(x%(x) + h.c 

where ^-fl( x ) = ^ / Tl( x )Xm(x)/v / 2 is the bare atom-molecule coupling. The molecule- 
molecule and atom-molecule interactions also follow from the above procedure but will 
be neglected in the following, since under the experimental conditions of interest 
the density of molecules is very small. 



3. Ladder summations 



For an the application of the hamiltonian in (|j) to realistic atomic gases we have to 
include all two-body processes, because at the relevant low densities three and more- 
body processes can in first instance be neglected. These are most conveniently included 
by a renormalization of the bare potential Vff(x — x') and the bare coupling g^i(x). 
Moreover, the molecules acquire a self energy. 

The interaction potential of the atoms in principle renormalizes to the many-body 
T(ransition) matrix. This renormalization is determined by a Bethe-Salpeter equation 
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Figure 1. Diagrammatic representation of the many-body T matrix. The solid 
lines correspond to single-atom propagators. The dashed lines corresponds to the 
interatomic interaction V^. 



which, within the for our purposes sufficiently accurate Hartree-Fock approximation 



171, reads 

TT 



T MB (k,k\K,z) = Ktfk-k') 



1 [l + iV(e K /2+k" - £t) + iV(e K /2-k" - fj) 
t7 2-~t 



V k „ Z — 6K/2+k" — CK/2-k" 

xT MB (k",k',K,z) , (5) 

where N[x) = \er x — is the Bose distribution function of the atoms, \x their chemical 
potential, and the inverse thermal energy. This equation describes the scattering 

of a pair of atoms from relative momentum k' to relative momentum k at energy z. Due 
to the fact that the scattering takes places in a medium the many-body T matrix also 
depends on the center-of-mass momentum K. The energy of a single atom is equal to 
ek = h 2 k 2 /2m. In the Hartree-Fock approximation the energy in principle also contains 
a mean-field correction which we neglect in first instance. We come back to this point 
in section [|, however. The diagrammatic representation of ([|) is given in figure p]. 

For temperatures not too close to the critical temperature we are allowed to neglect 
the many-body effects |17|], and (||) reduces to the Lippmann-Schwinger equation for 
the two-body T matrix. The effective interaction between the atoms thus becomes 
T 2B (k, k', z — €k/2). For the realistic conditions of the atomic gases under consideration 
here, i.e., small external momenta and energies, the two-body T matrix is independent 
of momentum and energy and equal to Ana^h 2 /m, with the s-wave scattering length 
of the potential Vj-f(x — x'). Therefore, we conclude that the renormalization of this 
potential is given by 

VWx-xO-^l^x-x'). (6) 
m 

The renormalization of the bare atom-molecule coupling is determined by the 
equation 

S MB (k, K, z) = g n (k) + W TMB (^ k '> K > *) 

v k' 
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Figure 2. Diagrammatic representation of renormalization of the bare atom-molecule 
coupling g n . 



1 + iV(e K /2+k' -fi) + N(e K/ 2-k> - n) 



x 1 ^uflO , (7) 

2 — CK/2+k' — CK/2-k' 

and is presented diagramatically in figure ||. Neglecting many-body effects, the coupling 
constant becomes g 2B (k, z — £k/2) with 



^ B (k,z)=^ n (k) + -5:r 2B (k,k' ) 

K k' 



(8) 



2! — 2ek' 

For the relevant small momenta and energies we are thus lead to the substitution 

<7u(x-x') -^(x-jO , (9) 

where g can be related to experimentally known parameters as follows. The resonance 
is characterized experimentally by a width AB and a position i?o- More precisely, the 
s-wave scattering length of the atoms as a function of magnetic field is given by 

(10) 



a(B) = a bg ( 1 



B-B Qt 

where ay, g denotes the so-called background scattering length. To make correspondence 
with the experiment we thus have that Off = a^ g . In order to reproduce the 



experimentally observed width of the resonance we have that g = h^2irab g ABAfi/m, 
since an elimination of the molecular field shows that ab g AB = mg 2 / \2tiTi 2 Aji) . 
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Figure 3. Diagrammatic representation of the self energy of the molecules. 



The self energy of the molecules, shown diagramatically in figure H, is given by 



^ MB (K^) = |£su(k) 



1 + iV(e K / 2 +k - At) + N(e K/2 - k - \i) 



T r / j 3 | 4, V — / 

v k Z — CK/2+k — CK/2-k 

xg MB (k,K,z). (11) 

In first instance we neglect again many-body effects which reduces the self energy in 
(0) to hJ: 2B (z-e K /2) with 

hZ 2B (z) = (Xm\V n Gv(z)V n \xm) , (12) 
where the propagator G^{z) is given by 

G n (z) = —L- , (13) 



z — 



with the hamiltonian 

^2 



H n = — + V n =H + V n . (14) 



We insert in a complete set of bound states \(f) K ) with energies E K and scattering 
states |0| c + ' ) ). The latter obey the Lippmann-Schwinger equation 

I0 k +) ) = |k) + — r^-V T M +) >, (15) 

where e k = e k + iO denotes the usual limiting procedure. This reduces the self energy 
KL 2B {z) to 

hT 2B {z) = £|<X m |Vn|0 K >| 2 ^- (16) 

K Z K 



(2r) JIU ™ " :-2e k ' 

where we replaced the sum over the momenta k by an integral. Using (P) and the 
Lippmann-Schwinger equation we have that 

^ B (k,2e+) = i=(x m |V r ul4 +) )- (17) 
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As a result we have for the retarded self energy hT>( + \hu), i.e., the self energy hT? B (z) 
evaluated at the physically relevant energy z = hu + that 

3/2 

- -g 2 ^iVh^ - (AfiB + E m ) , (18) 

where we have denoted the energy-independent shift, that results from the sum over 
bound states and the principal- value part of the integral in (|16|), in such a manner that 
the position of the resonance in the magnetic field is precisely at the experimentally 
observed magnetic field value B . 



4. Molecular binding energy and density of states 



Putting the results of the previous sections together, we find that the Bose gas near a 
Feshbach resonance is described by a coupled set of equations of motion for the atomic 
and molecular Heisenberg operators 



^ a (x,t) = exp 
^ m (x,t) = exp 



i 



-Ht 



ipl (x) exp 
^ m (x) exp 



4* 



(19) 



given by 



ih 



ih 



<9^ a (x,t) 

dt 



ft 2 v 2 

2m 



+ T b 2 s B ^(x,t)V; a (x,t) 



^a(x,t) 



+2#l(x,t)^ m (x,t) 

&l/} m (xL,t) 



dt 
-9 



h 2 v 2 



Am 



+ 5(B(t)) 



2 m 3 / 2 Id h 2 V 2 

— »t\hh— h 



(20) 



where the detuning is defined by 5(B) = Afi(B — B ) and = A-na^h 2 jm. This is the 
most important result of our work. Note that the time-derivative and gradient terms 
appear exactly such that both equations of motion are manifestly galilean invariant. 
Note also that an external trapping potential can just be added to the right-hand sides 
of these equations if required. It is interesting to mention that the above equations 
can immediately be generalized to a Fermi gas near a Feshbach resonance. Moreover, 
a simple Hartree-Fock approximation to the resulting theory reproduces exactly the 
interesting crossover physics recently discussed by Ohashi and Griffin on the basis of 
the Nozieres and Schmitt-Rink formalism for the normal phase of the gas |i~4| |. Having 
made this observation, it is now clear how the same crossover phenomena can be studied 
in the superfluid phase of the gas. Work in this direction is in progress and will be 
reported elsewhere. 

From d20|) we determine the retarded Green's function of the molecules 
G^(x, t; x', t'). For fixed detuning, the poles of its Fourier transform determine the 
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bound-state energy. This Fourier transform is given by 

h 



(21) 



e k /2 - 5(B) + (g 2 m 3 / 2 /2irti i )i<Jh(J - e k /2 

with ek/2 = fr 2 k 2 /4m the kinetic energy of a molecule. For positive detuning the 
propagator has a pole with nonzero imaginary part, which shows that the molecule has 
a finite lifetime in this case. In first approximation the energy of the molecule is e m (B) + 
e k /2 with e m (B) = 5(B) and its rate of decay equals T m (B) = (g 2 m 3/2 /AtiTl a )^J 5(B). 
For negative detuning the molecular propagator has a real pole at hu = e m + ek/2, where 
the bound state energy of the molecule is given by 



e m (B) = 5(B) 



4 S 

g m 
8n 2 h 6 



\ 



1 



16nW 
g 4 m 3 



5(B) - 1 



(22) 



Close to the resonance the bound-state energy is, using flKp, found to be equal to 
e m (B) ■■ 



m[a(B)] 2 ' ^ 
It is important to realize that this last equation, which is well-known to be true for 
a weakly bound state in a single-channel scattering problem [p2 |, has now thus been 
proven to be also valid for the case of a multi-channel Feshbach resonance. 

The physics of ( ETj ) is best understood by considering the molecular density of 
states, related to the retarded propagator by 

1 



■Im 



7T 



G^t^ (k, 



to 



(24) 



For negative detuning, it has two contributions as shown in figure |]. The first comes 
from the bound state, i.e., the (dressed) molecule, the second comes from the two-atom 
continuum. It is the latter part of the density of states that incorporates into our theory 
the rogue dissociation process put forward recently by Mackie et al. . More explicitly, 
the density of states is found to be equal to 

1 



Pu 



1 



-5(huj — ek/2 



+ -9(hu - e k /2) — 
vr [hui 



1+^3/2/(4^3^1) 

(g 2 m 3 / 2 /2TTh 3 )^/huj-e k /2 



,(25) 



e k /2 - 5(B)} 2 + (g*m 3 /4Tc 2 h G )(hLU - e k /2) 
and can be shown to obey the sum rule 

d(hu) p m (k, u) = 1 . (26) 

We thus conclude that the wave function renormalization factor of the molecules is given 
by Z = 1/[1 + g 2 m 3 l 2 1 (47r/i 3 y / |e m |)], which close to resonance is much smaller than one. 
Physically, this implies that in this case the wave function of the molecule is strongly 
affected by the interaction with the continuum in the incoming channel and contains 
only with an amplitude \fZ the wave function Xm( x ) of the bound state in the closed 
channel. 
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Figure 4. Molecular density of states. 



The molecular density of states changes when an atomic Bose-Einstein condensate 
is present, due to the mean-field interactions with the condensate. In particular, in 
this case we expect that the two-atom continuum part of the density of states is less 
important because of the mean-field barrier that the two colliding condensate atoms 
have to overcome. Mathematically this comes about since, to include the Hartree-Fock 
mean-field shift of the energy of the atoms, we have to replace in (|2lD ^Tiuj — ek/2 

by ^Tioj — 2hT l UF — ek/2, where frE HF denotes the Hartree-Fock self energy of the 
noncondensed atoms due to their interaction with the condensate. 

The mean-field shift of the thermal atoms leads to a change of the bound state 
energy as well. In equilibrium we estimate the magnitude of this shift by approximately 
calculating the self energy frE HF from 

2g 2 



HF 



2n 



rp2B I 
J-hn "■ 



(27) 



for a given condensate density no. In this manner we have calculated the mean-field shift 
of the bound-state energy as a function of the magnetic field for the Feshbach resonance 
at B ~ 154.9 G in the \ f = 2; m f = -2) state of 85 Rb. In figure | the results of this 
calculation are presented for two experimentally relevant condensate densities, namely 



n 



1.1 x 10 



13 



cm 



and n = 5.4 x 10 13 cm 3 



As expected, the shift of the 



bound-state energy is largest for the highest condensate density and decreases away from 
the resonance. Although ( |2"7| ) is only a first approximation to calculate the Hartree-Fock 
mean-field energy, figure || shows that the shift in the bound-state energy is significant. 
A more thorough calculation of the mean-field effect on the bound-state energy and 
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Figure 5. Mean- field shift of the bound-state energy as a function of the magnetic 
field. The solid line corresponds to an atomic condensate density of no = 1.1 x 10 13 
cm -3 and the dashed line corresponds to jiq = 5.4 x 10 13 cm~ 3 . The dotted 
line corresponds to the molecular binding energy in vacuum. The calculations are 
performed for the resonance at Bo ~ 154.9 G in the \ f = 2; to/ = —2) state of 85 Rb 

EE- 



also on the position of the resonance is postponed to future work. With respect to this 
remark, it should be noted that, since the retardation time of the interaction is only 
of the order of h/e m (B), the mean-field effects are not instantaneous. Therefore, the 
relevance of the mean-field shifts to nonequilibrium situations, such as in the recent 
experiments with time-dependent detuning [I], |2T| , is also subject to further research. 



5. Mean-field theory 



The mean-field theory for the coupled atomic and molecular condensates is found by 
taking the expectation value of (|20|) . We assume for simplicity that the temperatures 
are so low that the thermal cloud of the Bose gas can be neglected. It is, however, 
straightforward to include the mean-field effects of the thermal cloud in the same manner 
as in the by now standard Popov theory for weakly-interacting Bose gases. Furthermore, 
we neglect for simplicity the effect of the Hartree-Fock energy shift on the two-atom 
continuum and the molecular binding energy. The resulting mean-field equations are 
then given by 



ih 



ih 



dt 
dt 



^ s B |^a(t)|%(t)+2#:(t)^ m (t) 



6(B(t)) 




Vm(*)+#a(*) • 



(28) 



Note that in agreement with our remarks in the introduction the above equations 
contains no anomalous density or pairing field. One way to understand the reason 
for this is that the effects of the anomalous density are already included by using the 
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renormalized couplings and including the molecular self energy. Including these effects 
again would lead to double counting problems and, therefore, to a theory that is not 
gapless. Another way to understand it is that equation (|28| ) is explicitly U(l) invariant 
and the gaplessness of the theory is thus automatically guaranteed. 

A crucial ingredient in our formulation is the (nonlocal) term proportional to 
^Jihd/dt. It is this term that incorporates the correct binding energy of the molecules 
for negative detuning, and their lifetime if the detuning is positive. In principle this 
term must be treated as follows. Using the Green's function in (|2l|) we find that the 
wave function of the molecular condensate is, for time-independent detuning, given by 

= ? f dt'G^\t-t')^ a {t')+^{Q) , (29) 
n Jo 

where the Fourier transform of GW(0,w) is given by 



G(+\t-t') = -i6(t-t')Zexp 

iO{t -t')g 2 m 3 / 2 duo vT^e - ^* - *') 



-Umit-t') 

n 



« • (3°) 

■nh 6 Jo 2vr [hto - 5(B)} 2 + (g A nn? / ^ 2 h & )hu V ; 

This result can then be substituted in the equation for the atomic condensate wave 
function, which can now be easily solved numerically as we will show in future work. 



6. Conclusions and outlook 



We have derived from first principles an effective quantum field theory for Feshbach 
resonant interactions in atomic Bose gases. In future work we intend to apply this 
quantum field theory to study various equilibrium and nonequilibrium properties of 
ultracold atomic gases near a Feshbach resonance. This will include a further study of 
the mean-field shifts of the bound-state energy and the position of the resonance, as 
well as the study of the BEC/BCS crossover in the superfluid phase of a two-component 
Fermi gas |14]]- Moreover, we are now in a position to also study the normal phase for 
both for bosonic and fermionic gases. In order to apply the theory also to the recent 
pulse experiments with Bose-Einstein condensates of 85 Rb [p], ^] , we need to include a 
detuning that varies rapidly with time. This leads to some technical complications with 
the proper treatment of our nonlocal term that remain to be resolved. With respect to 
the latter experiments we also want to further study the importance of the quantum 



evaporation process, which previous work has shown not to be negligible [1 I 
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